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> ■ ABSTRACT 

£Nj ' We have examined X-ray and optical observations of two ultra-luminous X-ray sources, 

O X7 and X10 in NGC4559, using XMM-Newton, Chandra and HST. The UV/X-ray 

C*") 1 luminosity of X7 exceeds 2.1 x 10 40 erg s _1 in the XMM-Newton observation, and that 

of X10 is > 1.3 x 10 40 erg s _1 . X7 has both thermal and power-law spectral components, 
The characteristic temperature of the thermal component is 0.12 keV. The power-law 
components in the two sources both have slopes with photon index ~ 2.1. A timing 
analysis of X7 indicates a break frequency at 28 mHz in the power spectrum, while that 
• for X10 is consistent with an unbroken power law. The luminosity of the blackbody 

component in the X-ray spectrum of X7 and the nature of its time variability provides 
evidence that this object is an intermediate mass black hole accreting at sub-Eddington 
£5 ' rates, but other scenarios which require high advection efficiencies from a hollowed-out 

^ \ disk might be possible. The emission from X10 can be characterised by a single power- 

law. This source can be interpreted either as an intermediate mass black hole, or as a 
stellar-mass black hole with relativistically-beamed Comptonised emission. There are 
four optical counterparts in the error circle of X7. No counterparts are evident in the 
. error circle for X10. 

■ Key words: black hole physics — galaxies: individual (NGC4559) — X-rays: galaxies 

— X-rays: stars — accretion, accretion disks 



1 INTRODUCTION 

It has been known for some time that there are point-like, 
non-nuclear sources in some nearby galaxies, with luminosi- 
ties exceeding 10 39 erg s _1 (eg. Long, 1982, Fabbiano 1989, 
Colbert & Mushotzky 1999, Roberts & Warwick 2000, Mak- 
ishima et al. 2000) . The nature of these ultra-luminous com- 
pact X-ray sources (ULX) is still unclear. Fabbiano (1989) 
showed statistically that they are not foreground objects, or 
background AGNs. Some of the objects have been identi- 
fied with young supernova remnants (Fabbiano & Trinchieri 
1987), but many exhibit significant variability on long time 
scales. Others may be white dwarfs emitting via nuclear 
burning on their surfaces (Fabbiano et al. 2003). Low signal- 
to-noise spectra obtained by Chandra and ASCA show at 
least two classes of spectra: a pure power-law and a ther- 
mal blackbody-like spectrum (Strickland et al. 2001, Mak- 
ishima et al. 2000). They are most likely to be accreting 
compact objects (Colbert & Mushotzky 1999), but because 



their luminosity exceeds the Eddington luminosity of a 20 
Mq compact object (3 x 10 39 ergs s _1 ), sometimes by a fac- 
tor of ~ 15, this scenario suggests that the compact object 
may be an intermediate mass black hole (IMBH) of mass 
M.~ 100M Q . 

King et al. (2001) argued that the production of IMBH 
is difficult, especially when taking into account the con- 
straints on the mass-donor companion stars. However, they 
can be formed directly from the collapse of high mass stars in 
low metallicity environments (Heger et al 2003). ULX may 
also be numerous in elliptical galaxies (Angelini, Loewen- 
stein & Mushotzky 2001) although this is still unclear (Ir- 
win, Athey & Bregman 2003). An extreme case is their cre- 
ation as a primordial population in the first burst of star 
formation (Population III stars, Madau & Rees 2001). Al- 
ternatively, indirect channels may exist for the formation of 
IMBH, either through binary evolution (Fryer & Kalogera 
2001) or through the mergers of massive stars in Superstar 
Clusters (Portegies Zwart & McMillan 2002). IMBH may 
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also be formed through the mergers of lower mass black hole 
(BH) in cluster environments (Lee 1995, Miller & Hamilton 
2002). 

For higher metallicities, BH masses tend to be ^$20 Mq. 
Direct-collapse black holes (the most likely origin of 7-Ray 
Bursts) have probable progenitors which yield BH remnant 
masses with M.~ 14 - 23 M© (Fryer 1999). If other effects 
such as disk winds and magnetic fields are important, these 
masses could be lower. (For stars with insufficient angular 
momentum, weak 7-Ray Burst explosions or no explosions 
are produced, and the remnant can be much more massive, 
up to the mass of the star.) So, if ULX are powered by a stel- 
lar mass BH (which we define here as M, ^ 20M© ) , beaming 
or super-Eddington luminosities are required to explain the 
observed luminosities. 

Two different types of beaming can in principle ex- 
plain the luminosity of ULX in external galaxies. The beam- 
ing could be due to anisotropic emission of the radiation 
(King et al. 2001), or the relativistic motion of the radiation 
sources (for example Fabrika & Mescheryakov 2001). 

King et al. (2001) pointed out that if the accretion disk 
around the BH has a much lower scattering optical depth 
over a restricted range of solid angle (Afi) than in other di- 
rections, almost all the emitted X-rays will emerge in these 
directions. For a beaming factor b(= Af2/47r), the accret- 
ing BH could have a luminosity L — &L sp h, where L sp h is 
the luminosity inferred in the assumption of spherical accre- 
tion/isotropic emission. Such beaming is due to the special 
geometry of the accretion disk, thus it does not require the 
emitting plasma to be in relativistic motion relative to the 
observer. The spectrum of the emission could be thermal 
blackbody-like or a power-law. The chance of discovering 
such sources depends on the opening angle of the beaming 
cone, the duty cycle and the life span of the phases in which 
the mass transfer rate from the donor star is high enough to 
account for the observed X-ray luminosity. 

The other type of beaming is caused by the relativis- 
tic motion of the emitting plasmas. This would be the case, 
for example, if we were viewing the emission from BL Lac 
objects or the jets of SS 433-type systems at an angle al- 
most parallel to the jet orientation (Fabrika & Mescheryakov 
2001, Kording, Falcke, & Markoff 2002, Georganopoulos, 
Aharonian & Kirk 2002). The emission is Doppler-boosted 
by the relativistic motion, and it appears to be ultra bright, 
depending on the jet speed. As the emission is likely to 
be non-thermal direct synchrotron emission or (optically 
thin) Comptonised emission from relativistic beamed elec- 
trons, we would expect a power-law spectrum instead of a 
blackbody-like spectrum. The radiation should have strong 
linear polarization (even in the X-ray bands). The discov- 
ery probability of these type of sources is restricted by the 
opening angle of the jets. Fabrika & Mescheryakov (2001) 
estimated that about one such relativistic beamed source 
would be observed in each irregular or spiral galaxy, which 
is similar to the number of ULX per galaxy in recent obser- 
vations. 

It has proved difficult to distinguish whether ULX are 
stellar-mass Galactic BH binaries with beamed X-ray emis- 
sion or the more exotic IMBH systems. SS 433- and BL Lac- 
type relativistic beaming models can be excluded for sources 
with a thermal blackbody-like spectrum. Non-relativistic 
beaming (anisotropic emission) models will have difficulties 



if the X-ray observed luminosities are too high. Moreover, 
these kind of models would impose some signatures on the 
spectra either by Compton scattering or photoelectric ab- 
sorption. As the geometrical beaming model requires that 
the accretion disk is viewed almost face-on, any spectral 
features are expected to indicate a low inclination geome- 
try. If the geometrical beaming model is correct, the power 
spectrum of variability would be shifted to higher frequen- 
cies than in the IMBH model because the characteristic fre- 
quencies at the inner edge of the disk scale as 1/M. . For 
the relativistic beaming model, large amplitude variability 
should occur if the emission is direct synchrotron radiation. 

Further progress on ULX therefore depends on a de- 
tailed spectral study, a temporal study extending over as 
wide a range of frequencies as possible and a study of the 
location of the ULX in the host galaxies, in particular with 
respect to star forming regions. Ultimately, kinematic stud- 
ies of confirmed optical counterparts will prove to be defini- 
tive. 

It almost goes without saying that ULX are of signifi- 
cant interest. If they are IMBH, they are the link between 
stellar mass BH and supermassive BH. Their origin may 
be primordial, or the result of stellar or BH mergers. They 
will have a major impact on their surroundings. If they are 
beamed stellar mass BH, the beaming mechanism and the 
underlying structure of the disk are also astrophysically im- 
portant, and again, the scaling to AGN-generated jets is of 
interest. 



2 OBSERVATIONS 

This study concentrates on the study of two ULX in 
NGC4559, X7 and X10 (Vogler, Pietsch & Bertoldi 1997). 
Using the ROSAT PSPC, Vogler et al. (1997) identified 
seven X-ray sources falling within the D25 ellipse of the 
galaxy. The two brightest sources, X7, and X10, had intrin- 
sic luminosities ta 10 40 erg s" 1 in the 0.1-2.4 keV ROSAT 
band. X7 was located in the outer spiral arms, near an Hn 
region. They suggested that X7 could be a young supernova 
remnant in a dense cloud. X10 was thought to be associ- 
ated with the nucleus of the galaxy, although this source 
was noted to be extended, probably due to the superposi- 
tion of several sources, and offset from the optical nucleus 
by approximately 15 arcsec. 

Here we follow the Vogler et al. (1997) nomenclature; 
they are also known as XI and X4 respectively in Roberts & 
Warwick (2000), and 1X065 and IXO66 in Colbert & Ptak 
(2002). Roberts & Warwick (2000) commented that X10 is a 
nuclear source, while, from careful positional adjustments of 
ROSAT HRI observations, Pakull & Mirioni (2002) recorded 
it as a ULX displaced from the nucleus. We confirm here that 
X10 is non-nuclear. 

We make use of data obtained with XMM-Newton, 
Chandra and HST. The log of observations is given in Ta- 
blcQ The Chandra and HST data were retrieved from pub- 
lic archives. 

2.1 XMM-Newton data 

The XMM-Newton data were processed using SASv5.4. The 
XMM-Newton observations were centred on X7. The EPIC- 
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Figure 1. (a) The XMM-Newton EPIC-pn image of NGC4559. North is up and East is left. X7 is the bright source near the centre with 
the slightly fainter X10 located above and to the left. The inset shows the region around X10 and nucleus, with the Chandra positions 
superimposed, (b) The XMM-Newton OM UVW1 image of NGC4559. The locations of X7 and X10 are marked with arrows and small 
red circles and the centre of the galaxy is marked with a magenta cross. 



Observation 
Date 


Observatory 


Instrument 


Observation 
ID 




Duration 
(sec) 


Filters 


Mode 


1994-12-29 


HST (X10 field) 


WFPC2 


U29R2B01T- 


+02T 


160 


F606W 




2001-01-14 


Chandra 


ACIS-S 


2026 




9400 






2001-05-25 


HST (X7 field) 


WFPC2 


U6EH0505R- 


+08R 


2000 


F450W 








WFPC2 


U6EH0501R- 


+04R 


2000 


F555W 








WFPC2 


U6EH0509R- 


+0CR 


2000 


F814W 




2001-06-04 


Chandra 


ACIS-S 


2027 




10700 






2001-11-13 


HST (X10 field) 


WFPC2 


U6EA7101M- 


^02M 


460 


F450W 








WFPC2 


U6EA7103M- 


^04M 


460 


F814W 




2002-03-14 


Chandra 


ACIS-S 


2686 




3000 






2003-05-27 


XMM-Newton 


EPIC-pn 


0152170501 




40794 


medium 


full-frame 






EPIC-MOS1+2 






42130 


medium 


small-window (X7 only) 






OM 






41629 


UVW1 


image+(X7 only)fast 



Table 1. Log of the observations. 



MOS cameras were operated in small window mode to pro- 
vide sufficient time resolution, with the result that X10 was 
not seen in EPIC-MOS. The background was generally low 
during the observation, except at the beginning and for short 
episodes during the run. These intervals were excised from 
the data, reducing the effective exposure to 33.9 ksec for 
EPIC-pn and 36.8 ksec for EPIC-MOS. The image from the 
EPIC-pn camera is shown in Figure [Ua,. The 20 UVW1 (ef- 
fective band 2400-3200A) OM exposures each of 2 ksec du- 
ration were co-added (Figure^). 

An extraction circle of 30 arcsec radius was used for X7. 
For X10, we used a more complex extraction region of 20 arc- 
sec radius with an exclusion of 6 arcsec radius centred on the 
nucleus of the galaxy 13 arcsec away from X10. This reduced 
the contamination from the relatively fainter nucleus and 
another source 28 arcsec to the North and slightly West of 



X10 1 to a low level. Background was selected from source- 
free regions on the same CCD for the EPIC-pn data. For 
the EPIC-MOS X7 data, the small window mode prevented 
use of the same CCD for background subtraction, so regions 
from peripheral CCDs as close as possible to the source were 
selected. For the extraction, events were filtered according 
to event patterns 0-4 for EPIC-pn and 0-12 for EPIC-MOS 

1 In fact this source, which we call CXOU J123557.7+275807 is 
itself sufficiently bright to be classed as a ULX. With an absorbed 
power-law spectral fit, we estimate the emitted luminosity to be 
3 ± 1 X 10 39 erg s — 1 in the 0.3—12 keV band (approximately con- 
stant over the Chandra and XMM-Newton observations). Using 
an equally good diskbb fit, the inferred emitted luminosities are 
lower, ~ 1.5 — 2.0 X 10 39 erg s — 1 in the same band. Both this 
source and the nucleus are not resolved from X10 in Vogler et al. 
(1997). 
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(single and double events in each case) , and bad events were 
rejected using #XMMEA_PN=0 and XMMEA_EM=0. ar- 
fgen and rmfgen were used to create the response files. We 
checked that neither source is piled up. 

The extracted spectra for X7 from the three EPIC cam- 
eras are shown in Figure [5] We fitted the three datasets si- 
multaneously using XSPECvll.2 (Arnaud 1996). First we 
tried a simple power-law model, absorbed according to the 
tbvarabs (Tiibigen-Boulder variable absorption) absorber 
model with Solar abundance. This does not provide an ac- 
ceptable fit, xl — 1-8 Then we reduced the metal abundance 
to 0.1 of the Anders & Grevesse (1989) values. This pro- 
vided a relatively poor fit, xl = 1-12. Allowing the abun- 
dance to be fitted results in a lower xi = 398/381 = 1.04 
but with abundances < 0.01 of the Anders & Grevesse 
(1989) values, which would imply extremely low metallic- 
ities for the environment. We have no information on what 
the metallicities might be, but unless they are extremely 
low, this indicates that an additional soft component is re- 
quired. Adding a soft component improves the fit signifi- 
cantly, = 370/379 = 0.98, so that the component is sig- 
nificant at the > 99.99 percent level. Blackbody plus power- 
law (bb+po), disk blackbody plus power-law (diskbb+po) 
and Comptonised blackbody (bmc) 2 models, all with a tb- 
varabs absorber, all provide equivalently good fits to the 
data (in each case a WABS absorber model was also included 
to account for Galactic absorption). We found that the nor- 
malisation of the fits for the two EPIC-MOS cameras com- 
pared to the EPIC-pn is a factor ~ 0.93 of that expected, so 
we allowed this relative normalisation to be determined from 
the fits. The parameters derived from these fits are given in 
Table [5] There and elsewhere the units of normalisation in 
the fitted models are as defined in XSPEC. The low tem- 
perature of the thermal component means that we cannot 
distinguish between single and multi-temperature blackbody 
models in the EPIC band. The fit for the bb+po model is 
shown running through the data in the top panel of Figure^ 
while lower panel shows the x 2 residuals from the fits. 

The data from the three EPIC cameras were then co- 
added according to the prescription in Page, Davis & Salvi 
(2003). This spectrum with the fit from the parameters in 
Table [5] is shown in Figure [3] 

In the case of X10, an absorbed power-law model pro- 
vides an acceptable fit to the data and there is no need for 
additional components. The derived parameters are given in 
Table |3] and the spectrum is in Figure 01 

2.2 Chandra data 

In addition to the 2003 XMM-Newton observations, 
NGC4559 was observed three times by Chandra ACIS-S be- 
tween 2001 January and 2002 March (Table [T). On all occa- 
sions, the back-illuminated S3 chip was used. We obtained 
the Chandra dataset from the public archive and processed 
it with the CIAO software version 2.3. We inspected the 

2 The Comptonisation model BMC (Borozdin et al. 1999; Shrader 
& Titarchuk 1999) assumes a blackbody spectrum for the seed 
photon component, and is therefore preferable to the older 
COMPTT model which uses only the Wien tail of the blackbody 
spectrum. 
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Figure 2. The XMM-Newton EPIC-pn (blue, top), EPIC- 
MOSl(red) and EPIC-MOS2 (green) spectra of X7 simultane- 
ously fitted with the wabs(tbvarabs(bb+po)) model in Tabled 
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Figure 3. The XMM-Newton EPIC combined spectrum of X7 
(blue, top) fitted with the wabs(tbvarabs(bb+po)) model in Ta- 
bleland Chandra! (bottom, red) and Chandra2 (middle, green) 
spectra fitted individually with the wabs(tbvarabs(bb+po)) 
model in Table 151 

background count rates during the exposures, and chose to 
retain all three intervals in full. We applied the ACIS CTI 
correction thread, which also applied the newest available 
gain map. 

X7 and X10 are easily detected in all three exposures. 
Using the CIAO wavdetect task, we obtained image cen- 
troids and therefore X-ray positions. In addition, the nu- 
cleus of NGC4559 is evident in the first and third Chan- 
dra exposures, but is only at most marginally detected 
in the second. This corresponds to a reduction in emis- 
sion by a factor > 20 and a return to the previous level 
over periods of several months. We used the relative po- 
sitions of X10 and the nucleus to produce the extraction 
centroids for the XMM-Newton reduction above. We pro- 
vide the X-ray positions of X7, X10, the nucleus and the 
third brightest point source (CXOU J123557.7+275807) 
from this analysis in Table [I] The estimated uncertainty 
is 0.6 arcsec (90 percent confidence), derived from an ad- 
dition in quadrature of the centroid uncertainty, together 
with the systematic uncertainty of the spacecraft pointing 



ULX in NGC4559 5 



Parameter 



Value 



Parameter 



Value 



EPIC-pn, EPIC-MOS1, EPIC-MOS2 



EPIC-pn 



WABS(TBVARABS(BB+PO)) 


xl = 0.98 (369.6/379 dof) 


WABS(TBVARABS(PO)) 


xl = 0.85 (155.9/184 dof) 


njj (intrinsic) 


x 10 21 cm" 2 


njj (intrinsic) 


2.7tl'l x 10 21 cm" 2 


metal abundance 


0.311°;- 


metal abundance 


n QQ+0.53 
u - oo -0.19 


kT(blackbody) 


V.12+_IZ keV 


photon index 


2 05+ 07 


BB normalisation 


7.6±Jj x 10" 6 


PO normalisation 


1.97±°;^ x 10~ 4 


photon index 


9 90+O.O6 
z - zo -0.0S 


emitted flux (0.3-12keV) 


1.13 x 10~ 12 org cm -2 s _1 



PO normalisation 

MOS/pn normalisation 

emitted flux (0.3-12keV) 
blackbody flux (0.3-12keV) 
blackbody flux (0.01-12keV) 

WABS(TBVARABS(DISKBB+PO)) 

njj (intrinsic) 

metal abundance 

kT(inner) 

DISKBB normalisation 

photon index 

PO normalisation 

MOS/pn normalisation 

emitted flux (0.3-12keV) 
diskbb flux (0.3-12keV) 
diskbb flux (0.01-12keV) 

WABS(TBVARABS(BMC)) 

Tin (intrinsic) 

metal abundance 

kT(blackbody) 

photon index 

log(A) 

BMC normalisation 

MOS/pn normalisation 

emitted flux (0.3-12keV) 
flux (0.01-12keV) 



2.361°;?? x 10" 4 



,+0.05 
-0.03 



0.93Z u 

1.70 Xl0~ 12 ergcm~ 2 s -1 
0.43 Xl0~ 12 erg cm -2 s _1 
0.67 xlO -12 erg cm -2 s _1 

xl = 0.98 (372.7/379 dof) 
5.ll^3 x 10 21 cm- 2 



0.30 



3 

+0.27 
■0.20 
+0.006 



1 58+ 340 

9 90+O.O5 
z - z,3 -0.04 

2.3Tto.l° x 10' 



03+0 05 

U - JO -0.03 

1.97xl0~ 12 ergcm~ 2 s" 1 
0.72 X 10 -12 erg cm -2 s _1 
1.6 X 10~ 12 erg cm -2 s _1 



4.4+.Q7 X 10 21 cm" 2 



0.16 
0.12 
2.21 
0.06 



+0.27 
0.20 
+0.01 



+0.05 
0.05 
+0.08 
0.04 
+0.42 



1.041— x 10- 5 
03+0 05 



(see http://cxc.harvard.edu/cal/ASPECT/celmon/ ). This 
is confirmed by the coincidence between another Chandra 
point-source with a catalogued, isolated optical counterpart 
which has a clear counterpart in the USNO-B1.0 Catalog 
(Monet et al. 2003) within 0.4 arcsec (this source may be 
slightly extended). Another Chandra source is very close to 
the chip edge, so its astrometry may be affected, but even 
for this source the USNO-B1.0 Catalogue coincides to 0.16 
arcsec. A discussion of the X-ray sources other than X7 and 
XlO in this galaxy is left to further work. 

We used an extraction radius of 3.5 arcsec, with back- 
ground regions described by an annulus with inner and outer 
radii of 5 and 15 arcsec. We used psextract for extracting the 



Table 3. Parameters for the model fits to the XMM-Newton XlO 
data. The njj Galactic absorption for the WABS model is fixed at 
1.5 X 10 20 from Dickey & Lockman (1990). The metal abundance 
values are given relative to Anders &; Grevesse (1989) for the 
absorbing gas in the line of sight within NGC4559. All confidence 
intervals are 90 percent. 



Source 


RA 


Dec 


X7 


12 h 35 m 51!71 


+27°56'04'.'l 


XlO 


12 h 35 m 58?55 


+27°57'41'.'8 


nucleus 


12 h 35 m 57?64 


+27°57'35'.'8 


CXOU J123557.7+275807 


12 h 35 m 57!77 


+27°58'07'.'4 



Table 2. Parameters for the model fits to the XMM-Newton X7 
data. The njj Galactic absorption for the WABS model is fixed at 
1.5 X 10 20 from Dickey & Lockman (1990). The metal abundance 
values are given relative to Anders & Grevesse (1989) for the 
absorbing gas in the line of sight within NGC4559. All confidence 
intervals are 90 percent. 



Table 4. Positions for the three brightest point sources and the 
nuclear source in NGC4559 derived from the Chandra analysis. 
The estimated uncertainty is 0.6 arcsec (90 percent confidence. 



source and background, and building the response and aux- 
iliary response matrices. We also used the acisabs script pro- 
vided by the CXC to correct the auxiliary response function 
for the time-dependent decreased sensitivity of the ACIS-S 
detector at lower energies. 

Spectral fitting was possible for the first two Chandra 
observations. The third Chandra observation is too short 
for a spectral analysis, so we have used the countrates to 
provide only an additional point for the long-term intensity 
variations. 

The two background-subtracted Chandra spectra for X7 
are shown in Figure |3] below the combined XMM-Newton 
EPIC spectrum. Based on our experience with the 
XMM-Newton data, we have restricted our model fits to a 
BB+PO model, with the tbvarabs absorber model. The re- 
sults are given in Table |K| 

The two background-subtracted Chandra spectra for 
XlO are shown in Figureg|below the XMM-Newton EPIC-pn 
spectrum. As for the XMM-Newton data, we have restricted 
our model fits to a PO model, with a tbvarabs absorber 
model. The results are given in Table HJ and the three spec- 
tra and fits are shown in Figure 2] 

2.3 Variability 

2.3.1 Long-term variability 

From the fluxes derived from the best fit spectral parameters 
in Table|K|and|S| we obtained a long-term light curve for the 
0.3-10 keV band for both X7 and XlO. In the case of the 
third Chandra observation, we have derived the fluxes from 
the count rates assuming the average absorption, tempera- 
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Parameter 



Value 



Chandra 1 

WABS(TBVARABS(BB+PO)) 

njj (intrinsic) 

metal abundance 

kT(blackbody) 

BB normalisation 

photon index 

PO normalisation 

emitted flux (0.3-10keV) 
blackbody flux (0.3-10keV) 
blackbody flux (0.01-10keV) 

Chandra 2 

WABS(TBVARABS(BB+PO)) 

njj (intrinsic) 

metal abundance 

kT (blackbody) 

BB normalisation 

photon index 

PO normalisation 

emitted flux (0.3-10keV) 
blackbody flux (0.3-10keV) 
blackbody flux (0.01-10keV) 



xl = 1-16 (83.5/72 dof) 
3.S±li x 10 21 cm" 2 

0. 31 (from XMM-Newton) 

-,+0.06 
'-0.0 

8.5^;° x 10- 6 

1 sn+ - 08 

1. »u_ 08 

l.78+°; 2 ° x 10- 4 

1.66 x 10 -12 erg cm -2 s _1 
0.52 xlO^ 12 erg cm" 2 s" 1 
0.72 xlO -12 erg cm -2 s _1 



xl = 1-11 (105.8/95 dof) 

5.7±°;° x 10 21 cm" 2 

0.31 (from XMM-Newton) 

0.12±g;°i keV 

2.28lo5 x 10_5 
9 1O+0.08 
z - lo -0.08 

2.87+°; 2 ^ x 10- 4 



2.82x10" 



erg cm 
erg cm~ 



1.91 X 10 12 erg cm 2 : 



Table 5. Parameters for the model fits to the Chandra X7 data. 
The emitted fluxes are estimated in the 0.3—10 keV Chandra 
band. The njj Galactic absorption for the WABS model is fixed at 
1.5 X 10 20 from Dickey & Lockman (1990) and the metal abun- 
dance values are relative to Anders & Grevesse (1989) and fixed 
using the XMM-Newton value in Tablel2l All confidence intervals 
are 90 percent. 



Parameter 



Value 



Chandra 1 

WABS(TBVARABS(PO)) 

njj (intrinsic) 

metal abundance 

photon index 

PO normalisation 

emitted flux (0.3-10keV) 
Chandra 2 

WABS(TBVARABS(PO)) 
njj (intrinsic) 
metal abundance 
photon index 
PO normalisation 
emitted flux (0.3-10keV) 



5.6" 



0.86 (27.4/32 dof) 



-1.8 
-1.6 



0.33 (from XMM-Newton) 



1.99 



+0.22 
-0.21 



+0.27 
-0.21 

6.6 X 10 — 13 erg cm" 



L99 -o:2i x 10 ~ 4 



Xz< 
3.8 H 



0.33 (from XMM-Newton) 



1.08 (71.2/66 dof) 
;-° x 10 21 cm" 2 



1.82 



+0.12 
-0.13 
-0.26 
-0.23 
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Figure 4. The upper plot shows the XMM-Newton (top, blue), 
Chandral (bottom, red) and Chandra2 (middle, green) spectra 
for XlO together with the spectral fit and residuals using the TB- 
varabs(po) model in Tabled Note the residuals around 2.9 keV. 
These are more clearly evident in the lower plot, which shows the 
residuals of the spectrum of XlO from the model fit in the region 
around 2.9 keV for the Chandral, Chandra2 and XMM-Newton 
observations from top to bottom, respectively. 



ture and photon index parameters for the first two Chandra 
observations in Table|^]andE] This adds a fourth point. The 
light curve is shown in Figure^] 

We note that there is a significant change in the spec- 
tral index of the power-law component in X7 over the three 
observations, while the temperature of the thermal compo- 
nent remains constant. The spectral index of the power-law 
in XlO remains constant. 

The nuclear source is always too faint for spectral fit- 
ting in individual exposures. We have analysed the co-added 
spectrum, which can be fitted by an absorbed power-law 
with photon index ~ 2.2 and «h - 3 x 10 21 . We have used 
these average parameters to convert from count rates to 
0.3 - 10 keV fluxes. 



Table 6. Parameters for the model fits to the Chandra XlO data. 
Emitted flux is calculated in the 0.3-10 keV band. As for X7, the 
line of sight Galactic absorption is fixed at 1.5 X 10 20 from Dickey 
& Lockman (1990). All confidence intervals are 90 percent. 



2.3.2 Variability on shorter timescales 

We examined the variability of X7 and XlO within the 
XMM-Newton observations using Fourier spectroscopy. For 
XlO, EPIC-pn events within the energy range 0.3-9 keV (to 
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Figure 5. The long-term light curves of X7, X10 (left axis) and 
the nucleus (right axis). The ULX fluxes are determined using 
the tbvarabs(bb+po) and tbvarabs(po) fits respectively in the 
0.3-10 keV band. The nuclear light curves assume a constant 
power-law model with photon index 2.2 and nn= 3 X 10 21 . No 
XMM-Newton point is available for the nucleus because of the 
contamination from X10 (Figure 0. In terms of counts the nu- 
cleus is a factor of ~ 10 fainter than X10, which is similar to the 
level seen in the Chandra! and 3 observations. The error bars 
indicate 90 percent confidence ranges for the fluxes. 



maximise the signal-to-noise ratio) and within 30 arcsec of 
the source were selected, excluding circular regions around 
the two nearby sources with radii 10.5 and 8.2 arcsec. The 
events were binned into a light curve with a sampling period 
of 73.4 ms, the time resolution of the EPIC-pn detector. The 
presence of gaps in a light curve can contribute significant 
contaminating power over a wide range of frequencies, and 
therefore must be handled properly. We accounted for the 
gaps by subtracting the mean count rate (0.21 ct s _1 for 
X10) from the portions of the light curve with valid data, 
and assigning a rate of zero otherwise, which results in a 
zero-mean light curve. The light curve was zero-padded out 
to the next highest power of 2. The Fourier transform was 
then computed and the squared coefficients were rebinned. 
The resulting power spectrum was normalized to be the frac- 
tional mean squared variability per unit frequency interval. 
The result is shown in Figure [(J 

The spectrum appears to contain some variability at 
very low frequencies. To investigate further, we compared 
the results of fitting both a constant and a power-law model 
{P{v) oc v ~ a ) to the spectrum. The \ 2 f° r the power-law 
model was 10.5 for 11 degrees of freedom, compared to a \ 2 
of 25.9 with 12 degrees of freedom for the constant model. 
Using an F-test, we conclude that the variability is signif- 
icant at the 99.8 percent level. The best fitted power-law 
index was a — 0.9 ± 0.2, with a total r.m.s. variability of 
19 ± 7 percent. 

To see if overall background fluctuations could have con- 
tributed to the detected variability (such as low level flar- 
ing), we examined a 2800 arcsec 2 region of the same EPIC- 
pn CCD without bright sources. The mean background rate 
was 0.057 ct s _1 , or an expected background rate of 0.004 
ct s _1 in the X10 selection region. Since this rate is ~ 2 per- 
cent of the X10 mean rate, fluctuations of the background 
can be neglected, and the detected variability appears to be 
real. 

To investigate the variability in X7, we performed es- 



sentially the same procedure as for X10, with a few excep- 
tions. The EPIC-MOS1, -MOS2 and EPIC-pn events were 
selected within 30 arcsec of the source, and combined into 
one joint light curve, with a time sampling of 0.3 seconds 
(the time resolution of the Small Window image mode used 
for the EPIC-MOS observations). The mean counting rates 
were 0.29 ct s" 1 for the EPIC-pn, and 0.17 ct s" 1 for EPIC- 
MOS 1 and EPIC-MOS2 combined. The resulting power den- 
sity spectrum is also shown in Figure HJ We searched for 
apparent beats between the EPIC-pn and EPIC-MOS time 
sampling frequencies, but found none. 

Excess power is clearly detected, with an r.m.s. ampli- 
tude of 37 ± 2 percent. The improved detection significance 
is aided by the fact that the total count rate of X7 is almost 
a factor of two higher, and that the intrinsic variability itself 
is greater. There appears to be a break in the spectrum from 
flat to declining around 30 mHz. The precision is greatest at 
the break as a consequence of the binning and the change in 
slope (there are more Nyquist frequency bins per logarithmic 
bin as the frequency increases, so in general the statistical 
precision will increase as long as the signal is constant; how- 
ever, above the break, the signal decreases, so the statistical 
precision decreases). 

To treat this phenomenon more formally, we fitted both 
unbroken and broken power law models to the spectrum. 
We assumed a flat spectrum (qi = 0) below the break fre- 
quency, /{,. The unbroken power-law produced a \ 2 °f 81 
for 10 degrees of freedom, which improved to 14 for 9 de- 
grees of freedom when the broken power-law model was fit- 
ted, a significant improvement at the 99.99 percent level. 
Thus, we conclude that the break in the power spectrum 
is a real feature. The power-law index above the break was 
q 2 = 1.5 ±0.8. 

To test the flatness of the power density spectrum below 
the break, we performed an additional fit where the lower 
power-law index was allowed to vary. The best fit power-law 
index below the break was ot\ = 0.20 ± 0.14. The \ 2 value 
improves to 12.2, which is not significant enough to rule 
out an improvement purely by chance. Neither the upper 
power- law index, 02, nor the break frequency, ft, changed 
significantly with the new fit. 

Again, we investigated the fluctuations in the back- 
ground and found the background rate to be of order 2 per- 
cent of the source rate, and negligible. In addition, those 
background fluctuations that were marginally detectable 
were very low frequency (~ 10 - 5 Hz), very far from /{,. Fi- 
nally, we examined the power spectrum for the EPIC-pn 
and EPIC-MOS data individually, and found that similar 
variability was detectable in each, although at reduced sig- 
nificance. Thus, we are confident that the break is intrinsic 
to X7 and not the XMM-Newton detector systems. 

The best fitted break frequency was fb = 28 mHz. We 
determined confidence limits on fb by computing \ 2 on a 
grid of fb values, and allowing the other parameters to vary. 
The la confidence region for fb is 26-40 mHz. 



2.4 HST images 

The HST observations were retrieved from the pub- 
lic archive at ESO produced by the WFPC2 As- 
sociations Science Produces Pipeline (WASPP) (see 
http://archive.stsci.edu/hst/wfpc2/ ). Data for both the X7 
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Figure 6. The power density spectra for X7 (top) and X10 
(below) calculated from the XMM-Newton data. The best fitting 
model is shown with a shaded line. The Poisson noise level has 
been subtracted assuming that there is no intrinsic source signal 
above a frequency of 1 Hz. 



and X10 fields are available in three optical filters (Tabled . 
For X7 the images from each filter were combined to produce 
a single colour image, while for X10, we have used only the 
two later exposures. These images are shown in the upper 
panels of Figure 7. 

We found that the coordinate systems within the 
WFPC2 files were inconsistent from filter to filter by up to 2 
arcsec. We have therefore identified bright stars within the 
four WFPC2 CCDs corresponding to entries in the USNO 
A-2 and STSci Guide Star catalogs to recalculate the astro- 
metric solution. Only a small number of stars were available 
in each case, limiting the correction to a translation. Never- 
theless by this means the images can be registered to within 
0.3 arcsec. 

This allowed us to compare the Chandra and HST im- 
ages. Details of the images in the vicinity of the Chandra 
positions are shown in the lower left panels of Figure 7. Also 
shown is the error circle calculated by adding in quadrature 
the positional uncertainty from the HST images with that 
from the Chandra observations (which dominates) . This pro- 
cedure makes clear the environments of the two ULX. For 



X7, the 4 brightest stars in the error circle have 23 ^ V S; 
24 and 23 ^ R £ 24. There are no stars detected in the X10 
error circle with R ^$ 23.5 and B ^ 24.5. More details of this 
analysis will be presented in Soria et al. (in prep). 

The field of X7 was also observed at much lower spa- 
tial resolution by Roberts et al. (2002) using the Integral 
imager on the William Herschel Telescope. Their images are 
inverted (North is at the bottom) with respect to ours. We 
also note that there is a small displacement of ~ 0.7 arcsec 
between their position determination and our error circle. 
Mirioni (2003) shows Ha images of both the X7 and X10 
fields. 



3 DISCUSSION 

3.1 Long-term light curves 

Fabbiano et al. (2003) argued against an IMBH interpre- 
tation for the source CXOAnt J120151. 6-185231. 9, which 
has similar luminosity to X7, on the grounds that the wide 
variation in luminosity precluded it (assuming a constant 
temperature). Neither X7 nor X10 show a range more than 
a factor two in our data (Figure but our sampling is of 
course limited. Roberts et al. (2002) show the X7 light curve 
over a longer period including ROSAT points, together with 
some Chandra points. We have retrieved the ROSAT data 
from the archive, but find that the relative scaling between 
these and the Chandra and XMM-Newton points is uncer- 
tain by a factor of ~ 2 because of the smaller energy range in 
ROSAT, the lack of spectral information in the ROSAT HRI 
data and the confusion with CXOU J123557. 7+275807. We 
have therefore not incorporated them. Nevertheless, between 
Figure|S]and Roberts et al. (2002), the long-term variability 
of X7 is evidently modest. 

As noted in Section one possibility for the nature of 
ULX is that they are supernova remnants. Such variability 
as is observed in (Figure |SJ, and its non-monotonic nature, 
excludes the possibility that X7 and X10 are supernova rem- 
nants. 



3.2 Luminosity 

Assuming a distance of 9.69 Mpc (Sanders et al. 2003) and 
a Solar abundance, we obtain from the fluxes measured in 
the XMM-Newton data in (Table |SJ that the isotropic lumi- 
nosity of X7 in the 0.3-12keV band is 1.9 x 10 40 erg s _1 for 
the (bb+po) model. For X10, the 0.3-10 keV isotropic lumi- 
nosity (from Table[(jJ is 1.3 x 10 40 erg s . To estimate the 
total UV/X-ray luminosity of X7 we integrated the flux of 
the thermal blackbody component down to 10 eV. For the 
power-law component, we retained the 0.3—12 keV range 
from our observation, to avoid introducing arbitrary cutoffs. 
This gives a conservative lower limit to the contribution of 
the power-law component. The total UV/X-ray luminosity 
for X7 for the bb+po model is therefore 2.1 x 10 40 erg s _1 , 
of which ~ 8 x 10 39 erg s _1 is from the blackbody compo- 
nent. Using the diskbb+po model increases the luminosity 
to ~ 3.2 x 10 40 erg s" 1 , of which ~ 1.8 x 10 40 erg s" 1 is 
from the disk blackbody. For X10, we have only the power- 
law component, and so the conservative lower limit of the 
isotropic UV/X-ray luminosity is 1.3 x 10 40 erg s _1 . 
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Figure 7. HST images of the environments of X7 (3-colour) and X10 (2-colour) with insets showing the 90 percent error circles which 
are 0.7 arcsec radius. 



These luminosities are uncertain to the extent of the ac- 
curacy of the distance to NGC4559 and of the applicability 
of the spectral models. The Sanders et al. (2003) distance 
of 9.69 Mpc for NGC4559 is calculated from the redshift, 
using the cosmic attractor model in Mould et al. (2000). 
The value is the same as that calculated by Tully (1988). 
However, the distance to NGC4559 has not been measured 
with primary distance indicators. The intrinsic motion of the 
NGC4559 relative to the Hubble flow may affect the distance 
derived from redshifts, even when the correction has been 
made for our Galaxy. NGC4559 is a member of the Coma 
II cloud (Tully 1988), for which the only Cepheid distance 
measurement is that for NGC4725 at 12.7 Mpc (Gibson et 
al. 1999). The distance of 9.7 Mpc should therefore be con- 
sidered as a sensible lower limit. As noted above, the lumi- 
nosity/flux for X7 that we obtained depends on the assumed 
spectral models for the soft component. Nevertheless, all of 
the models yield luminosities exceeding 10 40 erg s _1 (in the 
0.3-10 keV band), and in some epochs the total UV/X-ray 
luminosity reaches ~ 4 x 10 40 erg s _1 . The derived total 
UV/X-ray luminosities of X10 are approximately half those 
of X7, from 8 x 10 39 erg s _1 (for the first Chandra observa- 
tion) to 1.5 x 10 40 erg s _1 (for the second Chandra observa- 
tion). X7 and X10 are therefore among the most luminous 
ULX known. These luminosities also require very high mass 



accretion rates, M ~ 1.7 X 10 _6 (0.1/»7) M yr" 1 , where 
rj{= L/Mc 2 ) is the efficiency of converting gravitational en- 
ergy into radiation. 

The X-ray luminosities of X7 and X10 clearly exceed 
the Eddington limit for 20 Mq compact objects with spher- 
ical accretion. This implies that both X7 and X10 could be 
IMBH if the emission is isotropic and the accretion is spher- 
ical. However, when these conditions are not satisfied, for 
example through inhomogeneity in the accreting material 
induced by instabilities such as those discussed by Begel- 
man (2002), it may not be necessary to invoke an IMBH to 
explain the observed luminosity. 



3.3 The nature of the compact accretor from the 
X-ray spectra 

3.3.1 The blackbody component 

The spectra of X7 show a soft thermal blackbody com- 
ponent, thus its bright X-ray luminosity cannot be easily 
explained by models assuming relativistic beaming of syn- 
chrotron radiation. The effective temperature of the black- 
body component is T e s ~ 0.12 keV. By assuming a "spheri- 
cal" geometry, a BB temperature of 0.12 keV and a luminos- 
ity of 0.8 x 10 40 erg s -1 (Table |5J, together with a distance 



10 



of 9.7 Mpc, imply a characteristic radius of 1.8 x 10 cm for 
the emission region. 

The nature of this characteristic radius is open to in- 
terpretation. If we assume that the thermal blackbody radi- 
ation originates from a compact opaque electron scattering 
sphere around a black hole, then the radius corresponds to 
the photospheric radius R p h of the scattering sphere. King 
& Pounds (2003) considered a wind/outflow model for such 
scattering region, in which the characteristic size R p ^ (which 
is defined by the photosphere), and the effective tempera- 
ture T e ff are expressed in terms of the mass outflow rate M out 
and the BH mass M, . By eliminating M out in their equa- 
tions (17) and (18), and inserting the values for R p h and T c s 
derived from our observation, we obtain M.~ 800Mq. 

We note that in this scenario, the emission is heavily ab- 
sorbed because of the dense wind material surrounding the 
system, as seen in Cyg X-3 (Heindl et al., 2003). Moreover, 
we also expect such wind material will be illuminated by the 
powerful X-rays from the accreting source, but no evident 
nebulosity is seen in the optical HST image (Figure 7) or in 
the narrow-band images in Mirioni (2003) and Roberts et 
al. (2002). 

Alternatively, the radius can be interpreted as the in- 
ner radius (R lB ) of the an accretion disk from which the 
soft thermal component originates. In this case, we may esti- 
mate the BH mass directly from the parameters of the fitted 
spectral models. The diskbb model, gives f?i n (cos i) 1 ^ 2 = 
1.21q'5 x 10 9 cm. If we assume that the accretion disk is 
truncated at the innermost stable circular orbit, then M,~ 
1A±1% x 10 3 (cos i)- 1/2 M Q for a Schwarzschild BH. From the 
BMC model, we obtain, instead, an effective radius R B ff and a 
colour tempertaure T co i (Schrader & Titarchuk 1999). By as- 
suming a hardening factor T co \/T e g = 2.6, a value applicable 
for BH systems, such as GRO J1655— 40, T e g together with 
R c a implies a BH mass ~ 1.7 ± 0.7 x lO 3 (cosi)- 1/2 M (see 
Schrader & Titarchuk 2003 for details). 

We note that all these arguments generally yield M, 
~ 10 3 Mq with L ~ 0.lL E dd- If we take these masses at face 
value, X7 is a strong IMBH candidate. 

The characteristic radius is not necessarily related to 
the last stable circular orbit around a BH. The soft black- 
body component may arise from a larger region. For an ac- 
cretion disk truncated at a radius much larger than the in- 
nermost stable circular orbit, the same fitted temperature 
implies a smaller M. but a larger L/L E dd ratio. For exam- 
ple, a system with M.w 200M Q , L ~ L E dd, and an accre- 
tion disk truncated at an inner radius R ~ 60R g (where 
R g — GM,/c 2 is the gravitational radius of a BH) would 
also have explained the low characteristic temperature de- 
duced from the spectral fits. However, in this scenario the 
emission from the region within the truncation radius must 
be suppressed, and it can be contrained observationally by 
the flux of the Comptonised (power-law) component. In or- 
der for the emission from this region not to exceed the value 
observed in the power-law component, a large fraction, > 90 
percent, of the radiation must be advected. A very large M 
is thus required to explain the observed X-ray luminosity. 
The situation becomes worse when smaller M, (and higher 
L/LEdd) are considered. For a system with M.~ 20Mq and 
an accretion disk truncated at an inner radius R ~ 600R g , 
then L ~ 10L Ed d- In this case, the advection is required to 
be extreme, with the accreting matter radiating only ~ 1 



percent of its gravitational energy. Whether this case is vi- 
able requires further theoretical investigation. 

M. may be also reduced if we allow the radiation to be 
anisotropic. A modest beaming factor b ~ 0.1, reduces M, 
by a factor of 10. Only very extreme values of b can reduce 
the black-hole mass from 10 3 to < 20M Q . However, Madau 
(1988) calculated beaming factors from thick disk structures 
to be < 5, and Strohmayer & Mushotsky (2003) argued that 
the QPOs seen in the M82 X-l source precluded substan- 
tial beaming. This poses difficulties for anistropic radiation 
(non-relativistic beaming) models, such as that suggested in 
King et al. (2001). 

Our conclusion above rests mainly on the high lumi- 
nosity and low effective temperature for the thermal (black- 
body) component obtained in the spectral fits for X7. We 
are aware that the temperatures obtained from the bb and 
diskbb models may not be very reliable parameters to de- 
termine the BH mass. Ebisawa st al. (2003) found that some 
ULX (e.g. "Source 1" in IC 342) can be explained more self- 
consistently within the conventional framework of accretion 
onto stellar-mass (~ 10 — 20 Mq) BH, when a slim-accretion 
disk model (see Wararai et al. 2000) is used in the spectral 
fitting instead of the models derived from standard optically 
thick geometrically thin disks (Shakura & Sunyaev 1973). 
ULX such as Source 1 in IC 342 tend to have high fitted 
temperatures (~ 1 keV), for the thermal blackbody compo- 
nent. However, X7 in NGC4559 is probably different from 
these sources since its thermal component has a much lower 
effective temperature. Similar cool temperatures are seen in 
thermal components of NGC1313 X-l and X-2 (Miller et al 
2003). 

3.3.2 The power-law components 

The XMM-Newton spectra of X7 and X10 clearly show a 
power-law component. This eliminates the possibility of nu- 
clear burning on the surface of a white dwarf, as investigated, 
for example, by Fabbiano et al. (2003) for the Antennae 
source CXOAnt J120151. 6-185231. 9. 

By analogy with the spectral classification of Galac- 
tic BH systems, X7 and X10 may be considered in a soft 
and hard spectral state respectively. The conventional ex- 
planation for the power-law component in the X-ray spec- 
tra of these systems is Compton emission from hot electron 
clouds (Sunyaev & Titarchuk 1980; Pozdnyzkov, Sobol and 
Sunyaev 1983; Zdziarski 1985). The accretion disk provides 
the seed photons, which are up-scattered to the keV energy 
range by the energetic electrons in a hot corona of an accre- 
tion disk or in a hot cloud surrounding the BH event horizon. 
The spectral index of the power-law depends mainly on the 
scattering optical depth and the temperature of the ener- 
getic electrons (for a thermal cloud). The photon index of 
the power-law component of X7 is T £s 2, corresponding to 
spectral index fsl,a value common to AGN and Galactic X- 
ray binaries. For modest scattering optical depths, say r w 1, 
this spectral index implies effective temperatures of about 
40 keV and 85 keV for the electron clouds with a slab and 
a spherical geometry respectively (Titarchuk & Lyubarskij 
1995). The electron temperatures will be higher for smaller 
optical depths. If we therefore extrapolate the luminosities 
of the power-law components to up to 50 keV, the total 
UV/X-ray luminosities for X7 and X10 rise to ~ 3.5 x 10 40 
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(for the (diskbb+po) fit) and ~ 1.7 x fO 40 (for the po) fit 
respectively. 

Alternatively, the disk photons can be up-scattered by 
non-thermal electrons, such as those in bulk relativistic mo- 
tion near the BH event horizon (Titarchuk, Mastichiadis & 
Kylafis 1996) or in relativistic jets (Georganopoulos et al. 
2002; Kording, Falcke & Markoff 2002). Photons Compton 
up-scatted by these non-thermal electrons also have power- 
law spectra. It can be shown that Comptonised emission 
from relativistic jets from stellar-mass BH can be Doppler- 
boosted and amplified by factors > 10 3 , making the lumi- 
nosity appearing super-Eddington (Georganopoulos et al. 
2002). However, such amplification can be observed only in 
a narrow viewing cone in which the jets are along the line 
of sight. Moreover, the observed luminosities are at the up- 
per limit of what this process can produce for systems with 
parameters similar to Cyg X-l. It may be suitable for X10, 
but not X7, with its high luminosity soft component. 

It has been suggested that the power-law X-ray com- 
ponent from BH binaries in the (low)-hard spectral state 
is due to direct synchrotron radiation (Fender 2001). If the 
observed X-ray spectrum from X10 is due to synchrotron 
radiation, the implied Lorentz factor 7 > 10 5 , and substan- 
tial power is also emitted in the optical and radio bands 
by an even larger number of electrons with lower 7. X10 is 
not listed in the FIRST radio catalogue (Becker, White & 
Helfand, 1995) (nor is X7), and its optical counterpart is not 
visible in the HST image of the field (Figure 7) at the level 
of R 23.5 and B ^ 24.5. Moreover, the large-amplitude 
variability on rapid timescales expected in the relativistic 
beaming synchrotron model for ULX (see Ulrich, Maraschi 
& Urry 1997) is not seen in the data of X10 (Figure [BJl. 
The power-law component is therefore unlikely to be direct 
synchrotron emission. 

In view of the similar properties of the power-laws of 
X10 and X7, their origins are probably the same, viz. they 
are both due to Comptonised emission from lower energy 
photons from the accretion disk. In sources such as X10, ev- 
idence of the soft thermal emission from the accretion disk 
can be hidden, if it emanates from radii somewhat larger 
than the innermost stable orbit and hence is softer. The soft 
disk emission is then overwhelmed by the harder Comp- 
tonised radiation from the hot electron clouds in the cen- 
tral region around the BH. If we assume that the accretion 
rate of X10 is sub-Eddington, the luminosity of the power- 
law component requires M, > 100 Mq. However, for super- 
Eddington accretion, a lower mass BH is permitted for the 
same observed luminosity. For M. < 20 Mq, the accretion 
rate must exceed 5 times the Eddington rate. 

To explore the maximum luminosity of a bb component 
consistent with the data, we included a thermal bb compo- 
nent with a fixed temperature of 0.12 keV . The best fit is 
obtained for no bb component, with a 90 percent confidence 
limit for the bb normalisation of < 1.76 x 10 -6 . The emitted 
flux from a thermal component in the 0.3-12 keV band is 
< 1.0 x 10~ 13 erg cm -2 s , < 8 percent of the total emitted 
flux in this band. The total emitted X-ray/UV flux in the 
bb component (0.01-12 keV) is < 1.5 x 10~ 13 erg cm~ 2 s~\ 
corresponding to a luminosity < 1.7 x 10 39 erg s _1 . Conse- 
quently, if there is a blackbody-like component in X10 with 
a temperature of 0.12 keV then it is an order of magnitude 
weaker than that in X7. Lower temperature thermal com- 



ponents would also be consistent with the data, as they fall 
outside of the XMM-Newton band. This would imply higher 
masses for the compact object, if the component is inter- 
preted as originating in an opaque scattering photosphere 
around it. 

3.3.3 Spectral residuals 

The spectral fit for X10 has residuals in the region of 2.9 
keV. These deviations are evident in the counts spectrum 
in Figure 2] Closer inspection shows that the residuals are 
in both Chandra and XMM-Newton spectra, but their na- 
ture is different: in the first Chandra observation there is a 
broad enhancement at this energy; by the time of the sec- 
ond Chandra observation, there is an absorption cut into the 
enhancement - the wings are still in emission, and by the 
time of the XMM-Newton observation there is only a narrow 
absorption feature. 

The origin of this feature is unclear. It may be instru- 
mental but we found no record of features at such energies 
in other Chandra and XMM-Newton observations. We there- 
fore proceed on the assumption that it is real. 

There are no strong spectral lines expected in this re- 
gion, the only candidates being S xvi at 2.62 keV and a 
fainter S xv line at 2.88 keV. Such lines are not expected a 
priori and are not observed in other bright BH candidates. 

We have considered whether X10 may be a background 
quasar. Supporting evidence would be the power-law spec- 
trum and the absence of any break in the power spectrum 
(unlike that seen in X7). In this case the feature at 2.9 
keV might be a gravitationally-redshifted Fe Ka line. With 
z — 1.2, the luminosity in the 2-10 keV range (using the 
ASCA band for comparison) is ~ 1.1 x 10 46 erg s _1 for 
the XMM-Newton observation, and ~ 6 x 10 45 erg s _1 for 
the first Chandra observation when the line is strongest. 
QSOs with this luminosity are bright but not uncommon at 
1 < z < 2. Variability by a factor of 2 is also normal. The 
fact that the emission line becomes weaker or disappears as 
the QSO becomes brighter is also common (Baldwin 1977). 
However, we believe this explanation is unlikely, since at this 
X-ray flux level the typical optical brightness of a QSO is 
R ~ 17 — 19, well above the upper limit of sources in the 
HST images. From the fitted values of tih (Tables |3] and |SJ 
we may expect an intrinisic extinction Ay ~ 2 magnitudes 
corresponding to an extinction ~ 4 mag at A ~ 2700 A (R 
band in the observer's frame). Even if we take this into ac- 
count we would still expect an optical counterpart brighter 
than R ~ 23 (Comastri, 2003), which is not seen in the HST 
images. In addition, the space density of background sources 
at this flux level is 0.2 per degree 2 (Rosati et al 2002) so the 
probability of being projected against the nuclear regions of 
NGC4559 is small, ~5x 10~ 4 . Finally, luminous QSOs do 
not generally exhibit X-ray emission or absorption features 
(Nandra et al. 2000). 

An alternative is that the feature is a Doppler-shifted 
line intrinsic to the source. If it is a red-shifted line, the only 
potential candidate is the Fe Ka lines, and then v/c = 0.66. 
Redshifted lines can originate from infall or a receding jet. 
However, either of these would have difficulty in produc- 
ing the narrow absorption features in the Chandra 2 and 
XMM-Newton data, when the central source is expected to 
dominate the continuum emission. If it is a blue-shifted line, 
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there are a number of possibilities of lines in the 1—2 keV 
range. Blue-shifting can be produced by outflow or an ap- 
proaching jet. It is puzzling that only one line is seen. This 
may imply an abnormal abundance of particular elements or 
a particular ionisation/radiation process. However, the tb- 
varabs fits to the XMM-Newton data do not indicate any 
abundance anomalies. The progressive change of the feature 
from emission to absorption, evident in Figure [I] could be 
explained in term of the jet/outflow scenario by condensa- 
tion and cooling of the outflowing material as its distance 
from the central X-ray source increases. 

In X10, the upper limit on the equivalent width for a 
narrow Fe K line at 6.4 keV is 80 eV, corresponding to a 
flux of 4 x 10" 7 ph cm s . The upper limit on a broader 
line depends on the fitted line model: for a gaussian line 
of width a = 0.5 keV, the upper limit on the equivalent 
width increases to 200 eV. In X7 the upper limits on the 
equivalent width for narrow and a — 0.5 keV Fe K lines 
at 6.4 keV is 130 eV and 45 eV respectively. It is not clear 
what the absence of an Fe K line means. In the Galactic 
BH the line strengths are variable and upper limits on the 
equivalent width of this order are frequently obtained. As 
an example, in the case of Cyg X-l, the upper limit on a 
narrow line from ASCA is ~ 30 eV (Ebisawa et al. 1996). 

3.4 The Nature of the compact component from 
the timing analysis 

Based on the similarity of the power spectrum in Galactic 
BH binaries and in AGN, many studies have investigated the 
possible link between the characteristics of the power spec- 
trum and the AGN mass (for recent examples see Markowitz 
et al. 2003, Uttley, McHardy & Papadakis 2002, Czerny et 
al. 2001). The power spectra are described by broken power- 
laws with characteristic break frequencies which appear to 
be related to the mass of the compact object. This is re- 
inforced by the extrapolation to the Galactic BH domain, 
where a scaled break frequency of ~Hz is compatible with 
M.~ IOM0, derived using binary kinematic arguments in 
Cyg X-l (Hasinger & Belloni 1990, Nowak et al. 1999a). 

The calibration of the M,—fb relationship is still rel- 
atively insecure. At lower masses, the relationship rests 
mainly on that for Cyg X-l (although breaks at ~Hz are 
also seen in other BH candidates) , while for the higher mass 
end, estimates of AGN masses are based on a relatively small 
sample determined mainly from reverberation studies. The 
physical connection between the two scales is also based on 
general timescale arguments (models include bright knots 
in the accretion disk - Lehto 1989, and Comptonisation - 
Kazanas, Hua & Titarchuk 1997), but more work is required 
to examine this in detail. For example, the cause for the 
observed range in power-law slopes is not clear, and other 
Galactic BH binaries such as GX339-4 show more com- 
plex power spectra, which have been modelled with QPOs 
(Nowak, Wilms & Dove 1999b). Nevertheless, for several 
AGN there are consistency checks (for example Edelson & 
Nandra 1999, Shih, Iwasawa & Fabian 2003), and the scaling 
of variability timescales is considered to be broadly valid. If 
this is accepted, then it is possible to use the power spectra 
in Figure |S] to infer information on the masses of X7 and 
X10. 

The situation is, however, complicated in that the power 



spectrum of Cyg X-l in the low-hard state is seen to have 
two breaks (Hasinger & Belloni 1990). One, fb2, occurs at 
a few Hz, where the slope flattens from a ~ —2 to ~ — 1, 
and a second, fbi, occurs at a 0.04-0.4 Hz (depending on the 
accretion rate, Hasinger & Belloni 1990) where the slope is 
consistent with being flat a ~ 0. In the soft-high state there 
appears to be no fbi (Reig, Papadakis & Kylafis 2002), and 
the slope continues as a ~ —1 to low frequencies. Some 
AGN have also been observed with double cutoffs (for ex- 
ample see Markowitz et al. 2003 and Uttley et al. 2002). 
The M. — fb relationship depends on which of the cutoffs is 
used at both high-and low-mass ends. Based on the index 
changes, Markowitz et al. (2003) suggest that the breaks in 
their sample correspond to the high frequency fbi breaks. 
In this case M.= 37// b M . (Markowitz et al. (2003) ex- 
cluded Akn 564 from their fits as they considered the break 
in this source at ~ 7.3 day to be a low frequency fbi break. 
However, it fits relatively well within their M.= 37/ /jMg 
scaling relationship, yielding a predicted mass of 2x 10 7 M Q , 
which is within a factor of ~ 2 of the reverberation mapping 
mass of 0.8 x 10 7 Mq.) If on the other hand, these are con- 
sidered to be the low frequency fbi breaks, as in the fits to 
the knee model in Uttley et al. (2002), then the relationship 
isM.= l// 6 M . 

We find for X7 (Figure © that f b = 0.028 Hz. The for- 
mal fits to the slope cti = 0.20 ±0.14 (Section ^. 3. 21 suggest 
that this cutoff is the lower frequency one fbi . On the other 
hand, the product P{fb) x fb ~ 0.02, which is comparable 
to the same product, ~ 0.01, measured by Markowitz et al. 
(2003), between the lower and upper break frequencies. It 
appears not to be possible to unambiguously select which 
cutoff to adopt. The indicated masses are M,= 38M© if 
fb = fbi an <l 13OOM0 if = fbi- If the break is instead to be 
interpreted as a QPO, then we note that the frequency is a 
factor ~ 10 lower than the "low frequency" QPOs seen in the 
stellar-mass BH candidates GX339-4 (Rieg et al. 2002) and 
RXTE J1118+480 (Belloni, Psaltis & van der Klis 2002). 

It is evident from Figure HJ that any break in the power 
spectra of X10 is at frequencies < 10~ 4 Hz. The scaling laws 
above would imply masses of M.~ 10 4 and 10 5 Mq in the 
two cases. 

3.5 The ULX environment 

Fabbiano & White (2003) discuss the environments of ULX 
in some detail. They note that in some cases ULX are asso- 
ciated with star forming regions, while in others, such as the 
Antennae, the ULX are offset from any star forming regions 
(Zezas & Fabbiano 2002). In some cases, such as NGC5204 
XI and M81 XI, an optical counterpart to the ULX has 
been identified (Roberts et al. 2002, Pakull & Mirioni 2002, 
for the former, Liu, Bregman & Seitzer 2002 for the latter): 
there the counterparts appear to be early-type stars. 

Outr HST images permit a detailed inspection of the 
environments of X7 and X10. For X7 there are five candi- 
dates within the error circle, of which the faintest and bluest 
is visible only in the F450W filter. One other similarly faint 
and blue candidate is on the edge of the error circle. Of 
course it is not necessarily the case that any of these stars 
is the true counterpart to the ULX secondary, since this 
may be, for example, a giant (which would not be visible) 
rather than supergiant. For X10 there are no obvious optical 
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counterparts in the error circle. However, this may be partly 
accounted for by the shorter exposure and coarser sampling 
of the X10 field in the HST observation. A detailed exam- 
ination of the environment of the two ULX and any effects 
of reprocessing of the X-ray luminosity will be presented in 
a companion paper (Soria et al, in prep). 

As noted in Section IMBH can be produced in low 
metallicity environments (Heger et al 2003), such as may 
be the case at least for X7, or as a primordial population 
in the first burst of star formation (Madau & Rees 2001). 
Volontieri, Haardt & Madau (2003) have examined the sub- 
sequent histories of these objects as they are incorporated 
in larger galaxies through hierarchical merging. Many will 
be incorporated in the central supermassive BH, but others 
remain in the halos or are ejected by three- or more body 
interactions (see Valtonen et al. 1994). They estimate that 
for a galaxy the size of the Milky Way (such as NGC4559), 
there will be ~ 10 wandering IMBH with Al. 150- 1000M©. 
The presence of two sources of the nature of X7 and X10 is 
therefore not surprising in this scenario. The question then 
is the expected frequency of capture of companion objects to 
fuel the accretion luminosity. These could be existing stars, 
or, more likely, new stars created by the IMBH triggering 
new star formation. This may be the reason for the asso- 
ciation of some IMBH with star forming regions (Fabbiano 
et al. 2001), although this association is not entirely clear 
(Fabbiano & White 2003). 

Another channel discussed in Section 0for the produc- 
tion of IMBH, this time non-primordial, may be through 
the mergers of stellar mass BH produced in star clusters or 
super-star clusters. The cluster will harden through interac- 
tions (Lee 1995), especially as it is perturbed by the envi- 
ronment, which may also be changing as a result of any hier- 
archical merging taking place. Such cluster remnants could 
also carry within them lower mass stars which may be the 
source of the material for accretion, when these are cap- 
tured. This may circumvent some of the difficulties raised 
by King et al (2001). The motion of these remnants through 
newly available gas when galaxy mergers take place may 
also provide the gravitational perturbation for new star for- 
mation regions. The ULX may appear displaced from the 
centre of the star forming region, depending on the time de- 
lay between the perturbation and the development of the 
star cluster. 



4 CONCLUSIONS 

In these investigations of the ULX in NGC4559 we have a 
useful conjunction of circumstances: observations from three 
great observatories of two sources at the upper end of the 
luminosity range of ULX. The X-ray spectra obtained us- 
ing XMM-Newton are among the highest signal/noise ratio 
achieved for such sources, with minimal nuclear contam- 
ination. These provide strong constraints on the spectral 
model fits. One of the objects, X7, shows both a blackbody 
and a power-law component in the X-rays, while the other, 
X10, has only a power-law component. The high count-rate 
permits an investigation of the temporal properties of the 
source to short timescales. The X-ray imaging from Chan- 
dra permits an accurate localisation of the sources within 



NGC4559, and HST provides deep imaging in the vicinity 
of the X-ray sources. 

The nature of ULX sources has been vigorously debated 
(Section^. Our data strongly constrain the interpretation of 
these particular sources. We can immediately exclude super- 
nova remnants, or nuclear burning on a white dwarf surface. 
We have interpreted our data in terms of accreting BH. 

In the case of X7 there are two lines of argument on 
which we have based our conclusions: the luminosity of the 
blackbody-type component, and the temporal properties of 
the source. In the first line of argument, the fluxes are de- 
rived directly from the spectral fits. The UV/X-ray luminosi- 
ties of this component can be inferred to be 0.8 — 1.8 x 10 40 
erg s~ . It is difficult to construct a scenario where signifi- 
cant beaming (b < 0.1) of this component can be achieved 
even if the emitting material itself is collimated. An IMBH 
accreting at L ~ O.lLEdd provides a simple explanation for 
the observed luminosities and blackbody temperatures, and 
at the same time a consistent explanation for the power-law 
component in terms of Comptonised emission. The masses 
derived from the disk fitting, or simple blackbody area mea- 
surements indicate that M.~ 2 x 1O 3 M0. This could be 
modified, however, by the (modest) beaming and inclina- 
tion factors. On the other hand, a model for a lower mass 
BH (Aim < 20A/q) fed from an accretion disk truncated at a 
large radius cannot be ruled out, if most of the energy within 
this radius is advected with high efficiency (~ 99 percent). 

In the second line of argument, the power spectrum ex- 
hibits a break at 0.028 mHz. If the relationship between 
break frequency and compact object mass can be inter- 
polated between stellar and supermassive BH, then this 
indicates X7 contains an IMBH. The masses determined 
from the temporal characteristics are M,~ 40Mq or M.~ 
1.3 x lO 3 Af0, depending on the calibration used. 

In the case of X10, the absence of the blackbody-type 
component prevents any inference of the mass from disk pa- 
rameters. The BH mass can be determined only from the 
luminosity: if this is sub-Eddington then we infer an IMBH 
for the compact object, while if we allow super-Eddington 
accretion or relativistic beaming it can be < 20M©. The 
power-law component has the same slope as that for X7, and 
we interpret this to arise by a similar mechanism, Compton 
scattering. Among the beaming models, Doppler-boosted 
Compton emission from a relativistic jet can explain the 
high observed fluxes without invoking an IMBH. The power- 
law is unlikely to be produced by direct beamed synchrotron 
emission as there is no evidence of sufficiently bright optical 
and radio counterparts, and the large amplitude temporal 
signatures expected are absent. 
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